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ABSTRACT 

We discuss the oscillation features of j3 Cephei, which is a magnetic star in which 
the magnetic axis seems to be oblique to the rotation axis. We interpret the observed 
equi-distant fine structure of the frequency spectrum as a manifestation of a magnetic 
perturbation of an eigenmode, which would be a radial mode in the absence of the 
magnetic field. Besides these frequency components, we interpret another peak in the 
frequency spectrum as an independent quadrupole mode. By this mode identification, 
we deduce the mass, the evolutionary stage, the rotational frequency, the magnetic field 
strength, and the geometrical configuration of (3 Cephei. 
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1. Introduction 

Beta Cephei (B1III, V = 3.2, V^sini = 25 km/s) 
is the prototype of a group of early-type pulsating 
stars. Since many of the (3 Cephei-type stars show 
multi-periodic pulsations, it is expected that a lot of 
information can be extracted from the pulsations of 
these stars. Recently, Aerts et al. (1994) and Telting, 
Aerts, & Mathias (1997) investigated the line-profile 
variability of (3 Cephei in detail, by means of an ex- 
tensive data set of high-resolution, high S /N spectra 
obtained at the Observatoire de Haute Provence dur- 
ing one month, and deduced that the variations show 
periodicity with at least five frequencies. The frequen- 
cies obtained by them from the intensity variation in 
the Si III A 4574 line are listed in table 1. Note that 
three frequencies among them (/ 3 , f±, and f$) are sep- 
arated from the main frequency f\ — 5.250 day -1 by 
—2/6 day -1 , —1/6 day -1 , and +1/6 day -1 , respec- 
tively. Hence, we regard that the four peaks (and the 
small peak at /i+2/6 day -1 ) are a part of a frequency 
quintuplet with equal spacing. The central frequency 
/i has been identified as the radial pulsation mode by 
Aerts et al. (1994) and by Telting et al. (1997). The 
second highest peak in the periodogram, / 2 seems to 
be an independent mode, and has been identified as 
an I = 2 mode by Telting et al. (1997). 

Table 1: Observed frequencies of spectral line variations. 
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Besides the pulsations in the luminosity and in the 
velocity field, a magnetic field has been reported at 
several occasions in the past. For a compilation of 
these results, see Veen (1993). No clear periodicity 
was found in these old measurements. Very recently 
however, Hcnrichs et al. (1999) re-measured the mag- 
netic field and found a clearly variable field with a pe- 
riod of 12 days and a semi-amplitude of 90±6G. The 
UV wind lines of this star also reveal a periodic varia- 
tion with both a 6 day and a 12 day periodic compo- 
nent (Henrichs et al. 1993), the 6 day period clearly 



being the dominant variation. From these results, it 
is concluded that the magnetic field of (3 Cephei is 
oblique to the rotation axis. 

In this paper, we try to understand the cause of 
the quintuplet and to extract information on the evo- 
lutionary stage and the geometrical configuration of 
(3 Cephei by using the oscillation data and the most 
recent magnetic data. A previous similar study was 
presented in Shibahashi & Aerts (1998) but they still 
used the older magnetic field measurements. 

2. Oscillations of a Rotating Magnetic Star 

Henrichs et al. (1999) find a variable magnetic field 
with a mean value of only ± 5 G. We assume here 
that the dominant component of the weak magnetic 
field is dipolar. Being influenced by such a magnetic 
field, the eigenmode, which would be a pure radial 
mode in the absence of a magnetic field, is deformed 
to have an axially symmetric quadrupolc component, 
whose symmetric axis coincides with the magnetic 
axis. Hence, the eigenfunction at the surface is char- 
acterized by means of a superposition of the spherical 
harmonic with £ = m = and that of I — 2 and 
Tn = with respect to the magnetic axis (Shibahashi 
1994, Shibahashi & Takata 1995); i.e., 

[Yo> B , 4>b) + aY°(e B ,ct> B )] exp(iwt), (1) 

where a is determined by the magnetic field of the star 
and the unperturbed eigenfunction of the mode. Since 
we are assuming that the magnetic axis is oblique 
to the rotation axis, the aspect angle of the pulsa- 
tion axis varies with the rotation of the star. There- 
fore, the contribution of the quadrupole component 
of the eigenfunction to the apparent intensity vari- 
ation changes with time and produces a quintuplet 
fine structure with an equal spacing of the rotation 
frequency in the power spectrum. Mathematically, a 
spherical harmonic expressed in the spherical coordi- 
nates with respect to the magnetic axis {9b,<Pb) is 
written in terms of (2£ + 1) spherical harmonics of 
the same degree I with respect to the spherical coor- 
dinates {Si,, 4>l)'' 

ytwbAb) = E E (-ir'Cw^Uw 

m '=-tm"=-t 

Y e m "(e L ,^ L )eM-im,'nt), (2) 

where {d^ m ,(p)} and {d^, m „(i)} are the matrices to 
transform the spherical harmonics expressed in terms 
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of the spherical coordinates with respect to the mag- 
netic axis, to those expressed in terms of the spherical 
coordinates with respect to the line-of-sight. The ex- 
plicit form of {d^m'iP)} f° r (■ — 2 cari be found in 
Kurtz et al. (1989). The angles 0,i G [0°,180°] are 
the angle between the magnetic axis and the rotation 
axis and that between the line-of-sight and the ro- 
tation axis, respectively. The observable variation is 
obtained by integrating equation ([!]) over the visible 
disc (8l = _[0, 7r/2] and 4>l — [0, 2tt\) after substituting 
equation (|2) into equation ([!]). Since the integral with 
respect to 4>l becomes zero unless m" — 0, the am- 
plitude of the component at lj + m'fl is proportional 



to 4ra'^) d m'oW (Shibahashi 1986). If an additional 
quadrupole component of the magnetic field is taken 
into account, the cigcnfunction is deformed to have 
I = 1, £ = 3, and £ = 4 components as well and then 
the expected power spectrum becomes a nonuplet (a 
9-fold multiplet) rather than a quintuplet (Shibahashi 
1994, Takata & Shibahashi 1994). The additional 
outer side-components are due to the £ = 3 and £ = 4 
components, and hence their amplitudes are expected 
to be much smaller than the central five components. 

3. Deduction of tan/3 and tani 

If the star is rotating with a period of 6 days, then 
the power spectrum of the eigenmode, which would 
be the radial mode in the absence of the magnetic 
field, is expected to reveal a quintuplet with an equal 
spacing of 1/6 day^ 1 . This can be explained by the 
observed quintuplet fine structure of fx, fs, f±, fe, 
and an additional very small peak in the power spec- 
trum (which we call fe). The relative ratios of the 
side-peaks' amplitudes to the central peak amplitude 
depend on the strength of the magnetic field. On the 
other hand, the relative ratios among the side-peak 
amplitudes depend on the angles (3 and i. By analyz- 
ing them, we can determine these angles. If we write 
the amplitude of the component at to + m'£l as A m >, 
then (A 2 + A_ 2 )/(A + -A-i) is given by (Shibahashi 
1986) 
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A 2 + A. 
A, +A. 
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= — |tan/3tani| . (3) 



Substitution of the observed amplitudes into the above 
equation leads to |tan/3tani| ~ 2.2, where we have 
used A 2 = A-2 since the presence of a frequency peak 
at Jq is not well established. The factor tan/3tani can 
be independently estimated from the variation in the 



magnetic field strength. Let r be the ratio between 
the apparent minimum and the maximum of the field 
strength, 

-Bobs .nun / ,\ 

— • (4) 



B, 



obs,max 



Then it is related with the geometrical configuration 
by (Stibbs 1950) 



tan [3 tan i 



1 



(5) 



Substitution of the observed values gathered by Hen- 
richs et al. (1999) together with their error, r G 
[—2.5,-0.6], leads to tan/3tani G [— 00, — 2.3] or 
[4.1,oo]. This result is seemingly in contradiction 
with the one derived from the amplitude ratio, but 
we will show in the Sect. [5] that the two estimates are 
compatible with each other and point to almost the 
same geometry. 

4. Identification of the Evolutionary Stage 

The second highest peak f 2 in the power spectrum 
has been identified by Telting et al. (1997) as a mode 
with 1 = 2 and m = +1. Telting et al. (1997) have 
chosen the rotation axis as the symmetry axis of the 
pulsation in their work. If the magnetic effects on the 
pulsation dominate over the effects due to the Coriolis 
force, as applies here, the symmetry axis of the eigen- 
function is the magnetic axis rather than the rotation 
axis. Following Telting et al.'s (1997) mode identi- 
fication, we assume that f 2 belongs to a quadrupole 
(I = 2) mode, but we assume m = with respect to 
the magnetic axis. As will be shown later, the fol- 
lowing conclusion about the evolutionary stage is the 
same even in the case of m ^ 0. Under the influence 
of the magnetic field, this mode is no longer described 
by a single spherical harmonic of £ = 2 and m = 0, 
and it is deformed to have components of some other 
£. However, the fine structure of this mode is expected 
to be difficult to detect, since the amplitude of the / 2 
component itself is small. 

Since we have two independent frequencies, /1 and 
f 2 , we can identify the evolutionary stage of the star. 
Figure [j] shows the frequency variation of a 1OM 
star. The left panel shows the radial mode frequency 
variation, while the right panel shows the case of £ = 2 
modes. As the star evolves from the zero age main 
sequency, its radius increases and hence the frequency 
decreases. The situation changes during the contrac- 
tion phase when hydrogen is exhausted in the stellar 
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center, and the frequency increases again in the hy- 
drogen shell-burning phase. The frequency variation 
of the I = 2 mode is more complicated because the fx- 
gradient zone around the convective core grows with 
evolution as a consequence of conversion of hydrogen 
into helium. The frequencies of the modes trapped 
there become higher. When the frequencies of two 
modes happen to be very close, they never degener- 
ate but repel each other. This is known as "avoided 
crossing" (cf. Unno et al. 1989). 




model number model number 



Fig. 1. — Frequency variation of radial modes (left panel) and 
1 = 2 modes (right panel) of a IOMq star with stellar evolu- 
tion. The vertical scale is in unit of day -1 , and the horizontal 
scale shows the model numbers (Model 1 is the zero-age main- 
sequence star, and Model 33 is the model at the turning point 
on the HR-diagram.). The solid lines indicate the observed 
frequencies /i (left panel) and /2 (right panel). 

By sweeping out evolutionary models of various 
masses, we search for stellar models whose radial and 
quadrupole mode eigenfrequencies coincide with the 
observed frequency /i and /2- Figure ^| shows the se- 
ries of these models on the HR diagram, calculated 
by using Pacyhski's (1970) program [Xq — 0.73 and 
Zq = 0.02). The two dashed lines show the models of 
which the quadrupole mode has the same frequency 
as /2, and the solid lines show the models whose ra- 
dial mode frequency coincides with /i . (3 Cephei must 
be on one of the crossing points of the solid lines and 
the dashed lines. The candidate is (i) a ~ 18M Q star 
at the middle of the hydrogen core-burning phase, or 
(ii) a ~ 12Mq star near the turning point, or (iii) a 
~ 9Mq star at the late stage of the hydrogen core- 
burning phase. 




X=0.73 
Z=0.0Z 



4.6 4.5 4.4 4.3 4.3 4.1 

log T,„ 

Fig. 2. — The HR-diagram of 8-20 M stars. The tick marks 
on each of the evolutionary tracks indicate Models 1, 10, 20, 
■ • • . The stars whose radial mode frequency coincides with the 
observed frequency /i are connected with the solid lines, and 
those whose quadrupole mode frequency matches with /a are 
shown by the dashed lines. The crossing points of the solid 
and dashed lines are the candidate locations for f3 Cephei on 
the HR-diagram. 

5. Deduction of i and (3 

The oscillation mode /i is identified as the radial 
fundamental mode in the case (i) or (iii) in the pre- 
vious section. On the other hand, it is identified as 
the radial first harmonic in the case (ii). First, let 
us consider the case (iii). The radius of the star for 
the case (iii) is estimated from the stellar evolution 
calculation as R ~ 6.5i?0. Since in our consideration 
2"7r/Sl = 6 days and T4sini ~ 25 km/s, this means 
i ~ 30° or 150°. Combining this with the first esti- 
mate of |tan/3tani| ~ 2.2 deduced from the power 
spectrum, we obtain (3 ~ 75° or 105°. The estimate 
tan (3 tan i ~ —6.14 deduced from the recent magnetic 
field measurements results in (3 ~ 95°. We obtain that 
the solutions from both estimates of tan (3 tan i point 
towards an almost equal geometry. The orientation of 
the magnetic axis can only be derived from the mea- 
surements of the magnetic field. We conclude that 
the angle between the rotation axis and the magnetic 
axis amounts to some 100° in (3 Cephei in the case of 
scenario (iii). 

If we take the case (ii), the radius of the star is 
larger than the case (iii) (R ~ 8.5i? Q ), and hence i 
becomes smaller (i ~ 20° or i ~ 160°) and (3 becomes 
close to 90° (f3 ~ 80° or 100° if | tan j3 tan i\ ~ 2.2 and 
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(3 ~ 93° if tan/3tanz ~ —6.1). Though the case (ii) 
cannot be ruled out, we think the case (iii) is more 
likely because the radial fundamental mode is more 
easily excited. The case (i) seems unlikely because of 
a high mass required. 

In the case of to ^ 0, the frequency of the quadrupolc 
mode is shifted by the Coriolis force by mCSlcos j3 
(Shibahashi & Takata 1993). Here C is determined 
by the equilibrium structure and the eigenfunction, 
and it is of the order of ~ 0.1 for the low order p- 
modes of I — 2. Since |cos/3| ~ 0.17, the frequency 
shift due to the Coriolis force is so small that we do 
not need to change the conclusion about the evolu- 
tionary stage discussed in the previous section. 

We have adopted the case (iii), and have calcu- 
lated the theoretically expected power spectra and 
compared them with the observations. Figure [| was 
calculated with an assumption of £? bs.max = 90 G, 
i = 30°, and f3 = 95° (the upper panel) and /? = 105° 
(the lower panel). The magnetic field was assumed 
to be mainly dipolar with a 10% contribution from a 
quadrupole component. (Note that the combination 
of i and (3 is reversible.) Figure || resembles the ob- 
served power spectrum, and it implies that our identi- 
fication of the pulsation modes, the evolutionary stage 
of the star, and the geometrical configuration are rea- 
sonable. 

6. Discussion 

The UV line variation of (3 Cephei implies that the 
rotation period is either 6 days or 12 days with the lat- 
ter being more likely as the UV line equivalent-width 
reveals a variation with alternating deep and less-deep 
minima (Hcnrichs et al. 1993). Based on this obser- 
vation, Henrichs et al. (1999) prefer an oblique rota- 
tor model with a dipolar magnetic field and with a 
rotation period of 12 days, the more so since this is 
the main period found in their 14 recent observations 
of the averaged value of the magnetic field over the 
stellar disk. Such a model can explain the UV data 
if an equator-on view is assumed. However, if the 
rotation period is 12 days and the magnetic field is 
a pure dipolc, then the power spectrum of the line- 
profile variations must show a quintuplet of an equal 
spacing of 1/12 day -1 for i ^ 90° and ^ 90°, or a 
triplet of an equal spacing of 1/6 day -1 for i = 90° 
or (3 = 90°; a quintuplet fine structure of an equal 
spacing of 1/6 day -1 is in that case unrealistic. 

In order to get an apparent quintuplet of the spac- 




Fig. 3. — Theoretically expected power spectrum. The model 
is a 9Mq star at the late stage of the hydrogen core-burning 
phase. The mode fi is identified as the radial fundamental 
mode. The angle between the rotation axis and the line of 
sight is assumed to be 30° and the magnetic axis is assumed 
to be inclined to the rotation axis by 95° (upper panel) and 
105° (lower panel). The magnetic field strength is assumed to 
be Bobs, max = 90 G, and the field is assumed to be a dipole + 
quadrupolc field, and the latter contribution is 10%. 



ing of 1/6 day -1 with the rotation period of 12 days 
from the oblique pulsator model, we have to assume 
that the magnetic field is almost entirely quadrupo- 
lar rather than dipolar and choose an appropriate ge- 
ometrical configuration to give only the five compo- 
nents among the nonuplet fine structure of the vis- 
ible amplitude. The former condition is necessary, 
because, otherwise, the eigenfunction would have an 
£ — 1 component, which would induce a pair of peaks 
separated from the central peak in the power spec- 
trum by 1/12 day -1 . In the case of a pure quadrupole 
magnetic field, the eigenfunction is characterized by 
means of a superposition of the spherical harmonic 
with I = to = and those of I = 2 and to = and 
£ = 4 and to = (Takata & Shibahashi 1994). The 
i — 4 component induces a nonuplet fine structure. 
But, in the case of i ~ 90° or (3 ~ 90°, both of the 
amplitudes at lu ± f2 and at u> ± 3£! happen to become 
much smaller than those at ui ± 2fi and to ± 40, and 
the fine structure looks like a quintuplet with an equal 
spacing of 1/6 day -1 (see figure |j). The combination 
of 2tt/H = 12 days and V c sin i ~ 25 km/s indeed leads 
to i ~ 90°, and one might consider that this would 
be favorable to explain the apparent quintuplet fine 
structure. However, in the case of a pure quadrupole 
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magnetic field, the observed magnetic field strength 
should vary as 



where 



i? bs oc Pi (cos 6) oc 3 cos 2 6 — 1 , (6) 



cos 9 = cos /3 cos i — sin/3sini cosfii. (7) 



Then, in the case of i ~ 90°, the observed magnetic 
field strength is expected to vary with a period of 6 
days rather than 12 days, and this is in contradiction 
with the observation (Henrichs et al. 1999). 



quadrupole field 
i=90' /3=50° 



U <J.± -J.Cj -J.O , <J.°± U.U 

Frequency (day ) 

Fig. 4. — The same as figure 3 but for the almost quadrupole 
magnetic field and i = 90° and /3 = 50°. The value of i has 
been chosen so that it is consistent with R ~ 6.5i?0, 2n/Cl = 
12 days and V e sin i = 25 km/s. The value of (3 has been chosen 
so that r ~ — 1.4 with an assumption of a quadrupole magnetic 
field and i = 90°. The dipole field contribution is assumed to 
be only 0.01%. 

In order to solve the controversy about the rotation 
period of the star, it is highly desirable that numerous 
new magnetic field measurements be performed over a 
much longer time base than achieved so far. We note 
that the older magnetic field measurements pointed 
towards very different values of the mean field, rang- 
ing from 70 G up to 800 G (Rudy & Kemp 1978, Veen 
1993). The new data obtained by Henrichs et al. 
(1999), however, are of much better quality. It would 
be extremely important to confirm the results by Hen- 
richs et al. (1999) and to achieve a better precision of 
the strength and the geometry of the magnetic field. 
This would allow a critical evaluation of the current 
theory of pulsations in hot magnetic stars. 
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